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Abstract 

We investigate the possibility that the CP violation due to the soft supersymme- 
try breaking terms in squark mixing can give significant contributions to the various 
7 related parameters in B decays, different from those of the Standard Model. We 
derive the new limits on 12) ll,lb.,rr and on {S23)ll,lr,rr from the recent data 
on D^~D^ oscillation as well as those on B^-Bg^ oscillation. We show that, to- 
gether with all the other constraints, the currents limits on these parameters still 
allow large contributions to the CP violating phases in B'^-BJ^ as well as D^-D^ 
oscillations which will modify some of the proposed measurements of 7 parameters 
in CP violating B decays. However, the current constraints already dictate that 
the one-loop squark mixing contributions to various B decay amplitudes cannot 
be competitive with that of the Standard Model (SM), at least for those B decay 
modes which are dominated the tree level amplitudes within the SM, and therefore 
they are not significant in contributing to CP asymmetries in the corresponding B 
decays. 



PACS numbers: 14.80. Cp, ll.30.Er, 12.60.Jv, 14.40.Nd 



1 Introduction 



With the B factories producing physics at full steam, it is important to investigate criti- 
cally the possibility of distinguishing standard Kobayashi-Maskawa (KM) |1| model from 
the other alternatives of CP violation. In particular, one needs to investigate if there 
is any non-KM mechanism that contributes to the measurement of CP violating phases 
which within the KM context, are labeled as a, /? and 7 (or 0i,2,3 in another popular con- 
vention in the literature). In fact, it is interesting to note that some of the CP violating 
asymmetries in B decays, which within the KM mechanism are considered identical, may 
correspond to different numerical asymmetries in an alternative theory of CP violation. 

One of the leading extensions of the Standard Model is the supersymmetric version of 
the theory. The addition of supersymmetric partners, as well as the necessity of supersym- 
metry breaking creates large variations of such models. In this paper we shall assume that 
the spectrum of the supersymmetric extension is minimal in the sense that no additional 
supermultiplet is introduced beyond the usual two copies of the Higgs doublet needed for 
fermion masses. In this theory, the new parameters are the coupling constants related to 
soft supersymmetry breakings. It is well known that these soft breaking parameters can 
give rise to new sources of CP violation. 

Among these soft breaking parameters, the ones most relevant to CP violation in K or 
B decays are the dim-2 soft squark mixing parameters (the matrices M? for left-handed 
squarks, and, M?, M'jj for right-handed squarks), as well as the dim-3 trilinear scalar 
Yukawa couplings (Y^, Yj^) which, after the spontaneous symmetry breaking, creates the 
mixings between left- and right-squarks (the matrices Y^{Hu) and Y^{Hd)^ here we follow 
the notation of Ref. 0). There are many discussions in the literature about whether the 
new supersymmetric contributions can give large enough e and e' [0, 0, ^. 

In Ref. 0, it was pointed out that even though the natural value of {5^2) lr in generic 
models may be of order 10~^, its contribution is big enough to saturate the experimental 
value for e'j^. One may wonder whether it is possible to use (^12) Li?, to saturate both ex 
and in kaon system[0, ^. However it was pointed out that in order to saturate both, 
the absolute value of (5i2)Lij has to be about 3 x 10~^ which is larger than its generic 
value. 

On the other hand, in Ref.[^, it was pointed out that if one takes into account the 
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isospin breaking effect of the supersymmetric AS* = 1 box diagrams, it is possible to 
account for e' even using the lm{6f2)LL with a mild fine-tuning. Note that, as emphasized 
in Ref.0, while 6lrs should be generically small due to its SU{2) breaking character, 
6ll and Srr do not have to be small. This opens up the possibility to use (^i2)ll alone 
to saturate both e and e'. 

In Ref.0, it was pointed out that while it is possible for {6f^)LR to contribute to (3 
parameter in B decays due to the B^-B^ oscillation, however its generic value is typically 
too small to account for the recent data from B factories, sin 2(3 = 0.59 ± 0.14 ± 0.05 (for 
Babarp) and sin 2/5 = 0.99 ± 0.14 ± 0.06 (for Belle|]). 

The purpose of this paper is to address the issue whether these new supersymmetric 
sources of CP violation can give rise to asymmetries which are usually associated with 7. 
We will consider only those B or Bg decay modes whose amplitudes are dominated by tree 
amplitudes within Standard Model (that is, the modes without "Penguin pollution"). We 
call these decays, non-Penguin type. We shall briefly comment on modes with Penguin 
contributions later on. 

We first derive constraints on squark mixing parameters based on the most recent 
data on B.-B^ and D^-D° oscillations. Within SUSY, AB = 1 box dia grams also arise. 
These can contribute even to modes that have only tree level contributions within the 
SM. Using the constraints obtained on the squark mixing parameters, we show that the 
contributions of AB = 1 SUSY box diagrams to the these (non-Penguin) B decays are 
negligible. However, we find that SUSY contribution can give rise to CP asymmetries 
through either the initial state B^-B^ oscillation or the final state D^-D^ oscillation. We 
compare these asymmetries with the predictions of KM model. 

2 New Limits on Squark Mixings 

In a comprehensive paper [|I^], Gabbiani et. al. work out various limits on the flavor 
changing couplings {5lj)LL,LR,RR^ where q can be M(up-type), (i(down-type) or £(lepton), 
and ij are generation indices. The 5's are dimensionless parameters defined as 

{5ij)AB = {m]j)ABlm\, (1) 
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where AB and ij stand for the chirahty and flavor respectively. For our purpose we will 
need only {S12) ll,lr,rr and ll,lr,rr- The earlier limits on these parameters, given 
in Ref. [0, are summarized in Tables 1 and 2. 

Recently there are new measurements on the D^-D^ oscillation as well as on the 
B^^-B^ oscillation. They can be translated into new limits on these 5 parameters. In 
particular, for D^-D^ oscillation, the new measurements give ArriD < 0.461 x 10^^*^MeV 
(in Ref . |lTl| ) . For Bg-Bg, the data have so far not been able to set a solid upper limit 
on the oscillation frequency, AMg, due to the error in the measurement of the higher 
frequency region |12|. The combined data from LEP and SLD give AMb, > 15 ps~^ at 
95 % CL. However, a hint of oscillation is observed (with large error) around Aitlb^ of 
17 ps^^ [|12|. The new measurements from LEP and CDF can also be combined to give 
Ar,/r, = 0.16l[5:[]^, or Ar,/r, < 0.31, at 95% C.L. This can be combined with the lattice 



calculation of AVjAMs, = S.5Z\%t x lO-^jg in SM, to give Atub, = (291^?) ps"^ . 
In presence of SUSY contributions, we expect, this lattice estimate to change. To obtain 
reasonable limits on the ^23 parameter, we use the suggestive values of Attib^ = 8, 17,45 
ps~^ as typical value in our study. Note that these are not yet serious experimental limits, 
however as commented later, our physical conclusions on B decays in the next section, 
are not significantly altered even if AMb^ turns out to be one order of magnitude larger. 



X 












^\^m,,mnR\ 


0.3 
1.0 
4.0 


4.7 X 10-2 

1.0 X 10-^ 
2.4 X 10-^ 


6.3 X 10-2 
3.1 X 10-2 
3.5 X 10-2 


1.6 X 10-2 

1.7 X 10-2 
2.5 X 10-2 



Table 1: Limits on Re (5i2)yiB (^12)00 froni Amo, with A, B,C,D = (L, R), for an average 
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squark mass rriq = SOOGeV and for different values of x = m'^/m'^. 



X 


V ^-^J LL 


V ^-^/LR 


0.3 


4.4 


1.3 X 10-2 


1.0 


8.2 


1.6 X 10-2 


4.0 


26 


3.0 X 10-2 



Table 2: Limits on the 



"23 



from 6 ^ S7 decay for an average squark mass rUq — 500 GeV 



and for different values of a; = m?/m?. For different values of rriq, the limits can be 
obtained multiplying the ones in the table by (mq;/500 GeV)^. 
For -B°-S° oscillation, the {5l^)ab 

mixing contributes to the operator 

Qs = ^'fPl^rPI > 

Qb = slhisihl, (2) 

plus the operators Qi,2,s obtained from the Qi,2,s by the exchange L R. Here qr^l — 
|(1 ±75)g, and a and P are color indices. The color matrices normahzation is Tr(t^t^) = 
^S"^^- The AB = 2 effective Hamiltonian reads: 



^eff = -^_[{sls)^^{24Q,xMx) + 66Q,Mx 

+ {sl,)l^{2AQixMx) + 66QiMx)) 

+ K)^^ (4)^^ (504 g4 a; /six) - 72Qj,ix) 

+ 2AQ^xU{x) + 120 Qje{x)) 

+ {6l,)l^{204Q2xU{x)-36QsxU{x)) 

+ {6l,)l^ (204 g2 a: /6(x) - 36 ga feix)) 

+ (4)^^ (4)^^ (-132g4/6(x) - 180g5/6(:r 



(3) 
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where x = m'^/m'^, rriq is the average squark mass involved in the box diagram, rrig is the 
gluino mass and the functions fei^x) and fe^x) are given by : 

6(1 + 3x) \nx + x^ - _ g^, _^ 17 



6(x-lY 



~ , , 6x(l + x) In X — — Ox^ + 9x + 1 

Note that j^ix: = 1) = 1/20 while j^ix = 1) = —1/30, therefore they cancel a lot in the 
combination 24x/6(x) + 66/6(x) = —1 for x = 1. 

The matrix elements of the operators Qi, Q2 are defined as 

(5.|gi|i?.) = (5) 



5 .0 ,.0 Mi ,^ 

S-, 
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/I Ml^^s . (6) 



12' 

The matrix elements of the other operators in Eq.(|^) can be obtained in terms of those 
given in Eqs.(^J^). The bag factors B and Bs parameterize the non perturbative contribu- 
tions to the matrix elements and have been evaluated |T^ on the lattice to be i3 = 0.9±0.1 
and Bs = 1.25 ±0.1. 

Similar equations apply to the calculation of the mass difference of D^-D^ system. 



Using the same naive estimate of hadronic matrix element used in Ref.[]TU[, the results 
may be summarized in Table 3. Note that, while for the Bg, the hadronic matrix elements 
discussed above, have been obtained with some rigor, those for the D meson assume an 
universal bag parameter of unity. 

As one can see by comparing Table 3 with Table 1 and Table 2, the recent (and 
coming) data do provide significant improvement on the limits on (5i2)yiB- 
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X 



0.3 
1.0 
4.0 



X 



0.3 
1.0 
4.0 



12) LL 



2.58 X 10-2 
5.46 X 10-2 
1.28 X 10°^[^-i 



^ {sis) 



2 

LL 



0.16,0.23,0.38 
0.34,0.50,0.81 



0.80,1.17,1.90 



12)LR 



3.43 X 10-2 
1.72 X 10-2 
1.90 X 10-2 



5? (52^3) 



2 

LR 



0.20,0.29,0.46 
0.11,0.17,0.27 
0.13,0.18,0.30 



RR\ 



8.52 X 10-3 
9.49 X 10-3 
1.34 X 10-2 



'23 



LL 



'23 



RR 



0.06,0.08,0.13 

0.06,0.09,0.15 
0.09,0.13,0.21 



Table 3: Values of Re (^«i)cD' 'with A,B,C,D = {L,R). The upper part of the 

table are derived from saturating ArriD < 0.479 x lO-^'' MeV by squark mixing contribu- 
tion. The lower part of the table are based on suggestive values Am^^ = 8, 17,45 ps"^. 
We use an average squark mass rriq = 500 GeV and choose different values of a; = m'~/m'~. 
The constraints on (%)rr are the same as those on (Sij)^^ in the Table. 



3 SUSY contributions to 7 

The clean measurement of 7 has been a challenge, leading to several attempts at providing 
feasible techniques to measure it. SUSY contributions to 7 therefore vary, depending on 

the method used. We therefore first discuss the various methods proposed to measure 

(— ) 

7. The original suggestion [|15i for cleanly measuring 7 involved the decays 



and D^^pK^, where D^p stands for the CP eigenstate of D. However, since it is 

virtually impossible to tag the flavor of the D meson, the method was improved. Ref. [|l^ 

(— ) 

considered the produced, to subsequently decay to at least two final states. The mode 

(~) 

D^*K* was proposed in Ref. |T^. For the purpose of this paper the arguments 

(— ) 

made to the generic mode B"^ -^D^ applies to all the methods in Ref.|l5|, |16|, p!7| . 



Alternative modes involving Bg mesons have also been suggested to measure 7. They 
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include B^/B^^ DfK^ |[T|] (see Fig. 1) and its final state vector meson analogue 
B^/B'^ — s> D*^i^*^|l9|. One may note that all the modes discussed above only have tree 
level contributions. Other methods involving B Ktt together with B KK, etc., 
which include penguin contribution, have also been considered|^. However, they involve 
theoretical assumptions like the inherent use of SU(3) or factorization assumption. 

All CP asymmetries arise from a relative phase between two decay channels to the 
same final state. This may arise either due to two or more contributions to the direct 
decay or, due to the oscillation of the initial or final state neutral meson. 

Let us consider the CP asymmetry appearing in charged decay modes B^ — D^K^ 
or D'^K^. In SM, the first stage of the decay subprocess at quark level is either due 
to 6 — s> u{cs), which is severely Cabbibo suppressed but complex in the Wolfenstein (or 
Chau-Keung) convent ion pT|, or, due to 6 — > c{us) which is doubly Cabbibo suppressed 
and real in the same convention. By observing decays of and to a common final 
state / (which may be a CP eigenstate), one achieves the required interference. The 
relevant effective Hamiltonian for the quark level process b u{cs), in SM is 

l{gj/M^)V:^V^,iuLi,bL)isLi''cL) , (7) 

where V*gVub ~ AX^e~'''^] while for the quark level process b — > c{us), it is, 

ligl/M^)V:McLl,bL)i-SLYuL) , (8) 

where V*gVcb ~ AX^. The relative phase gives rise to CP violating parameter 7 directly. 
In the scenario of squark mixing, the effective Hamiltonian for b —* u{cs) is 

KV"^'^^ = -T^(^i2)^^(^23)ll {2^xUx) + 66h{x)) {u,^,bL){s,YcL) + ■ ■ ■ , (9) 

with X = m^/mi. We only list the contribution from the channel LL in chirality. Other 
amplitudes due to the insertion of other 5 parameters are easily obtained from Eq. (^. 

For SUSY contribution to the AB = 1 amplitude to be relevant to CP asymmetry, 
it has to be a sizable contribution to the decay amplitude. To estimate the SUSY con- 
tribution we take x = 1, 24x/6(x) + 66 f^i^x) = —1 and optimistically use ~ 0.02 
iKfel ~ 0.003, {S'^2)ll ~ 0.06, (6^3) LL ~ 0.7. The amplitude ratio of the SUSY to the SM 
is about 



(a2/108)10-^ GeV~^(100 GeV/MgY x 0.06 x 0.7 /lOO GeV^ 



{AGf/V2) X 0.003 V ^< 



X 10-2 . (10) 
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Assuming that the SUSY AB = 2 box diagrams dominate in the Bg — 5^ and — 
oscillations, the same ratio can be expressed more directly in terms of measured quantities 
as, 

fDfBsVM^/ V2 ^ ' ^ ' 

The numerical result is based on the central values of the parameters, AMd = 0.07 ps^^, 
fo = 0.2 GeV, and /^^ = 0.23 GeV, as well as the suggestive value AMb^ = 45 ps~^ If 
one takes Mg = 500 GeV in Eq. (|10D, The number comes up to be 10~"^ as in Eq. (p^T]) . 
Therefore even after taking the parameters to be favorable to the SUSY contribution, 
such one loop contribution is still much smaller than the highly KM suppressed SM tree 
amplitude. (Even if the color suppression of the contribution of the KM operator to 
B^ D^K^ is taken into account, SUSY is still very much a minor contribution to 
decay amplitudes). Therefore we conclude that SUSY contribution to the AB = 1 box 
diagram (Fig. 2) is irrelevant to the CP asymmetry as long as we limit our consideration 
to 5| and 6^2- The same conclusion can be applied to all the B decay modes. 

One may hence conclude that the only possible SUSY contributions to the CP asym- 
metries arise from initial state AB = 2 or final state AD = 2 transition, or both. Among 

the various methods to determine 7, B^ DfK"^ get contributions from initial state 

_^ (— ) _^ 

B^ — Bg oscillation, B"^ get contributions from final state — D oscillation 

(— ) _^ _o 

and and Bg </> get contributions from both B^ — B^ and — D oscillations. 

We first consider the contributions from initial state 5^ — 5^ oscillation. Since SUSY 

particles are heavy, we assume that is not modified by SUSY and parameterize the 

SUSY contributions by M^^^^ = M^^ye''^. Hence, we have, 

^ 12 ^ -1-12 I^^h2 ^ /^r,N 

s and are SM parameters, with s ~ 0(10^^) and 0^0. In terms of the SUSY 
parameters discussed earlier, for the LL chirality we have. 



2 

LL 



The expression for M^*^ may be taken from Ref . . Using 



AM = -2Re(^(Mi2-^ri2)) , (13) 
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the expressions for AM including SUSY may be written as, 



1 

2 , . \ 2 



AM = 2|Mf2*^| Re|^l + 2ycosr7 + ?/2 - ^ -is(cos0 + ycos(0-?7) 

= 2|Mf2^^| ^{l + 2y cosr^ + y^) (1 + 0{s^)) (14) 



The effect of SUSY on the meson mixing parameter {q/p) can be expressed as, 



(15) 



\pI Mfi' {l + ye^^) ■ 
Since, Mf2^ is real for 5^, the argument Ob^ of {q/p) is given by, 

es = - tan-^ ( -y'M^v)-^y-^-v ] . (^g) 

2 Y 1 + y"^ cos(2r7) + 2?/ cosr^y 

?/ may be solved using Eq. (|1^ as a function of r], leading to the oscillation phase as a 
function of rj. Using the upper limit of AMb^ = 45 ps^^ and theoretical estimate of M 



SM 
12 ' 



from Ref. |2^, we find that values of Ob^ are allowed in the full range (—45°, +45°). If, 
however, M^^^^ > Mf^^ i-e. in the limit, ?/ oo, we get Ob, = —i] = — Arg {{S23)aa)^ 
with A = L or R. 

The mode Bg DJK^ has been proposed to measure sin^7 (see Fig. 1). In the 
presence of SUSY contributions to 5° — 5^ oscillation, the angle 7 gets modified, 7^7' = 



^KM ± Gbs^ where, Ob^ is given by Eq.(|T^) and ^km is the contribution from Kobayashi- 
Maskawa phase. The contribution from SUSY to both the modes Bg —>■ DfK^ and 
Bg {il!/J)(j) are identical. This is very much different from the KM predictions in which 
the asymmetry in B^ ^ ip/ J + (p is negligible while that is 5^ DjK~^ is large. Note 
that the phase Ob^ can be measured directly using the mode Bg {ip/J)(f), which has a 
large branching ratio and should be easier to measure. 

In decays of the type B^ — > D^K^, as discussed earher, 7 is measured by utilizing a 
possibility of interference between the two quark level processes h —>■ cus and h —>■ ucs. 
Along the decay chain the c or c produce in the final state a or mesons respectively. 
The two contributions are added and interfere if both or decay to the same final 
state f£) and have a relative phase 7. Here, fo is one of the states that both D and D 
can decay into, such as K~ti'^ or CP eigenstates K^K~, vr+vr", Kgir^ or Kg(j). In fact 
the whole discussion can be applied to the modes in which final state is replaced 
by TT^ or In the presence of D^-D^ oscillation there are additional contributions to 
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this process. As discussed in details in Ref.[^, |2J], there are many different types of CP 
violation that can manifest themselves in such decays. SUSY can have potentially large 
contribution to the D^-D^ oscillation, even providing a large phase to the oscillation. 
One can estimate the D^-D^ oscillation phase, by repeating the procedure used to 
determine Ob^i except that here Vf^ /M^^ cannot be ignored. In the SM model, 7 is 
large, and 6d is small. However, if SUSY were to dominate, 6d = — Arg((5"2)yiB(5i2)cD), 
with A,B,C,D=(L,R), could give large contributions to CP asymmetries even if 7 is small, 
especially Arg(((5"2)iL or Arg((5]'2)/jH- Ref. has considered the various possible CP 
violating asymmetries in detail for arbitrary 9o- First of all, there is the phase in the 
decay amplitude, which in KM model, is exactly the 7 parameter. Then there is the phase 
in the D^-D^ oscillation, Od = aig^qn/pD) (here, qa pn are the composition amplitudes 
in the D-D system). In addition, the differences in strong final state phase shifts in 
decays (A^) and in D decays (A^)) may become relevant for some CP observables. 
There is CP violation of the type proportional to sin 7 sin A^ which is due to the CP 
asymmetry in the decays B"^ DK^, the final state phase shift is needed to produce the 
CP asymmetry for charged B decay as expected. There is CP violation proportional to 
sin 6'/) sin A/) which is similar to the CP asymmetry in the decays D fo- There is CP 
violation proportional to sin^/jcosA/j which is due to D-D oscillation. Finally, there is 
CP violation of the type proportional to sin (7 + Oo) cos A^ which is due to interference 
between B ^ D decays and the subsequent D-D mixing. Last category is of course 
most interesting. In the KM model, 7 is large, however 6'£) is small, while in the SUSY 
model we consider, 7 is small but 6*/) = aig{{5'^2)\A) '^^'^ t)e quite large (here, (^^4) can 
be either LL or RR). Of course, in the fourth type of CP violation listed above, Od in 
SUSY model can duplicate the effect of 7 in KM model, however as discussed in Ref . 



in details, there are enough CP asymmetries that one can measure to distinguish the two 
contributions in principle. 

Finally, we have so far avoided discussing B decays that can receive significant con- 
tributions from penguin or chromo-dipole moment types of diagrams such as, at quark 
level, h —>■ dss, b —>■ sss, b —>■ ddd or 6 ^ sdd, either in KM or in SUSY models. They 
contribute partially to i? — > tttt, 5 ^ A'tt, B Kef) and other processes. While the 
Ai? = 1 box diagram SUSY contribution is still negligible for these processes as long 
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as one considers only (Jfg and 6^2 parameters, as discuss in Ref . [p^ , the SUSY penguin 
and/or chromo-dipole moment types of contributions can play important role in their CP 



asymmetries. This issue is discussed in Ref. p6 



4 Conclusion 

We list below the results of our study of (5^^ rr)23 and (^^l _r/?)i2 on the CP asymmetries 
in Bd or decays. 

(I) We use the stringent D^-D^ mixing data to obtain tighter limits on {S^l rr)!^- 
Based on suggestive values of we illustrate the amount of improvement 

can be made on the constraint of (5^^ rr)23 in the near future. 

(II) For CP asymmetry in B decays, SUSY can give large contribution to Bs decays due 
to Bs-Bs oscillation, but cannot give large contribution to the complex phase in the 

decay amplitude via AB = 1 box diagram. Thus SUSY produces CP asymmetries 

(— ) 

different from those in the KM model. For example, while in KM the B^ 
process has large phase 7 and Bg {ilj/J)(j) has negligible CP violating phase, 
however in SUSY, both processes receive the same CP phase due to the Bg-Bg 
oscillation. 

(Ill) For Bd decay modes or charged B^ decays modes, the CP asymmetries due to the 

(— ) 

above 5's are negligible. However an exception is B^ -^D^ K^, where there can be 
SUSY contribution to the CP asymmetry due to the final state D^-D^ oscillation. 
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Standard Model Contribution 




Fig. 1 Quark diagrams for the tree-amplitudes in SM for the process (Bs/Bs) K^D 

Squark Mixing Contribution 




Fig. 2 The Box-graph contribution in SUSY via squark mixing to the process [Bg/Bs) 
K+D-. 
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